Scalar field collapse and cosmic censorship 
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We analyze here the final fate of complete gravitational collapse of a massless scalar field within general 
relativity. A class of dynamical solutions with initial data close to the Friedmann-Lemaitre-Robertson-Walker 
(FLRW) collapse model is explicitly given and the Einstein equations are integrated numerically in a neighbor- 
hood of the center. We show that the initial data space is evenly divided between the dynamical evolutions that 
terminate in a black hole final state and those that produce a locally naked singularity. We comment on the 
genericity aspects of the collapse end-states and the connection to cosmic censorship conjecture is pointed out. 

PACS numbers: 04.20.Dw,04.20.Jb,04.70 Bw 

Keywords: Gravitational collapse, black holes, naked singularities, scalar fields 



(N 

o 

oo 



O. 

cr 

i ■ 

$H - 



> 
oo 

(N 

\o 

O 

(N 



X 



The Cosmic Censorship Conjecture (CCC) 12] has been one 
of the most prominent open problems of modern general rel- 
ativity. The original idea of Penrose was that the evolution 
of a generic and physically realistic matter cloud, subject to 
Einstein equations, would never lead to spacetime singular- 
ities that are visible to faraway observers. Over the last few 
decades much work has been done on CCC, with multiple for- 
mulations being proposed towards a possible proof, and on 
what is meant by genericity of a certain collapse outcome, a 
concept for which there is no precise definition in gravitation 
theory today. 

While a proof of the CCC has not been found despite many 
efforts, and its very formulation has been questioned, many 
models of physically realistic gravitational collapse leading to 
the formation of naked singularities have been obtained and 
studied (see for example [2] and references therein). These 
models mainly deal with fluids such as dust, perfect fluids, and 
other more general matter fields, which have been used widely 
in astrophysical considerations. These are typically described 
by averaged properties where the microscopic structure of the 
fluid itself may have been neglected in the very final stages of 
collapse (see, however, JH]). In any case, it is not clear whether 
the fluid approximation will continue to hold till the ultra-high 
densities that characterize the singularity are reached, or what 
is the final form of matter and equation of state valid at such 
ultra-high densities. 

On the other hand, much work has also been done on scalar 
field collapse (H-Hl). 

Scalar fields are fundamental fields 
satisfying the Klein-Gordon equation and therefore, despite 
the fact that they have never been observed in nature, consti- 
tute a perfect candidate for a fundamental field whose validity 
will hold at any scale. Also, as shown by Choptuik jjj, the 
black hole threshold in the space of initial data for classes of 
massless scalar fields show universality and power-law scal- 
ing of the black hole mass, which correspond to critical phe- 
nomena. These are explained by the existence of exact so- 
lutions which are attractors within the black hole threshold 
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and are typically self-similar. The critical phenomena evolve 
a smooth initial data to arbitrarily large curvatures visible 
from infinity, via the formation of infinite time naked singu- 
larities that form at the boundary between the continual col- 
lapse to a black hole and the region where collapse reverses, 
turning to a dispersal (see [9] for details). The genericity 
and stability aspects of naked singularities arising in certain 
classes of collapsing massless scalar fields were investigated 
by Christodoulou [5], showing these to be non-generic within 
the context of a certain parent space. 

We investigate here complete collapse of a class of spheri- 
cally symmetric massless scalar fields within general relativ- 
ity. We show explicitly that there exist sets of non-zero mea- 
sure of initial data, which are arbitrarily close to the FLRW 
homogeneous collapse scenario when the collapse begins, and 
which terminate in a naked singularity. We note that the naked 
singularities considered here are of a different nature com- 
pared to the infinite time singularity in the Choptuik like anal- 
ysis in the fact that in our models the spacetime is no longer 
assumed to be self-similar, and these are much similar to the 
central singularities of fluid models that may be locally visi- 
ble. 

The energy-momentum tensor for a massless scalar field 
4>(r, t) is given by, 
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We consider here the class of scalar fields which admit one 
timelike and three spacelike eigenvectors (see e.g. IToll ). and 
consider the continual gravitational collapse of such a field. 
The FLRW scalar field collapse is an example of such a sce- 
nario. Then T a b can be written in comoving coordinates and 
becomes diagonal. The spherically symmetric line element 
then depends on three functions of t and r (which are the 
comoving time and radius respectively), which are <?oo = 



2v(r,t) 



9n = e 2 ^ r '*) and 322 = #33 sin 2 
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As we can always choose here a frame along the eigenvec- 
tors of the energy-momentum tensor of the scalar field, we 
have To = 0. We then get cjxj)' = 0. We shall con- 
sider here the class of scalar fields for which <f>(t) ^ 0]. 
Then the components of the energy momentum tensor become 
T h a = ^e~ 2v 4> 2 [— 1, 1, 1, 1], from which we see that the mat- 
ter behaves like a stiff perfect fluid with an equation of state 
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p = p, with p given by the equation above. 

The system of Einstein equations in units where 8irG = 
c = 1, can be written as 



the definition of p given by the energy-momentum tensor we 



F' = pR 2 R' , 
F = -pR 2 R , 
= d t {R 2 e^-^ 
GR' = Iv'RG . 



(2) 
(3) 
(4) 
(5) 



where primed and dotted quantities denote partial derivatives 
with respect to r and t respectively, and we have defined 
G = e~ 2 ^'R' 2 . Here F(r, t) is the Misner-Sharp mass of the 
system describing the amount of matter enclosed within the 
shell labeled by r at the time t. 
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Note that equation together with the definition of the scalar 
field implies that <fi must obey the Klein-Gordon equation. 
The system of Einstein equations together with the equation 
of state and the Misner-Sharp mass is a set of six equations 
in the six unknowns p, p, ip, v, R and F, and is thus a fully 
closed system for this case. 

There is a scaling degree of freedom that we can use to set 
R(r, ti) — r. We then write R(r, t) = rv(r, t) and rewrite all 
the equations in terms of the variable v. Thus v(r, t) is such 
that v — 1 corresponds to the initial time and v = to the time 
of occurrence of the singularity, that is, the shell labeled by r 
becomes singular at the time t s (r) given by v(r,t s (r)) = 
( lITTll . lfl2ll ). The condition for continual collapse is then given 
by v < 0, while the situation where v changes sign at a finite 
time corresponds to a halting and dispersal of the scalar field. 
We can rewrite all the relevant functions in terms of the new 
coordinates {r, v}, therefore for any function X(r, t) we have 
X 1 = X, r + X, v v' where v'(r, t) itself is treated as a function 
of r and v. 

Avoidance of shell-crossing singularities, which are gen- 
erally believed to be a breakdown of the coordinate system 
rather than true singularities, is obtained if we impose R 1 > 0, 
which corresponds to v + rv' > 0. Since v is always positive, 
there is always a small neighborhood of the center r = for 
which no shell-crossing occurs, and our analysis holds good 
for such a collapsing cloud. 

In order to describe physically valid scenarios the energy 
conditions must be fulfilled throughout collapse. In this case 
these are always satisfied once p is taken to be positive, which, 
provided R' > 0, corresponds to F' > 0. Furthermore regu- 
larity must be imposed on the initial data in order for p and 
the metric functions to be well behaved. We shall require 
that the Misner-Sharp mass behaves like F — r 3 M(r, v) with 
M'(0, v) — in order for the matter density to be regular and 
without cusps at the center at all epochs of collapse earlier 
than the singularity. 

An important point here is, to solve the system of Ein- 
stein equations we write all equations in terms of the mass 
function M(r, v) and its derivatives. From equation (|4]i we 
get R 2 e^- U <j) = f(r) where f(r) is an arbitrary function 
related to the kinetic energy of the collapsing cloud. From 
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■v cj) /(2M_ V ), while using the equation of 
state and the above definition for / allows us to write G as 
G(r, v) = -v 2 {vM tV - 3M - rAf r ) 2 /2/ 2 Af „. Now using 
equation (f2]i and (fJJ we can write the differential equation for 



v 1 = W(r, v) 



vM, 



3M + rM, 



2Af , 



while from the Misner equation © we obtain 



V{r,v) 
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with negative sign taken in order to describe collapse. From 
the above we see that for collapse to occur we must require a 
'reality condition', namely + > 0. If this condition 
is not satisfied throughout the evolution then the system will 
reach v — in a finite time and collapse will then reverse to 
dispersal. Finally the whole system reduces to a second order 
partial differential equation in M and its derivatives, which 
can be written as 



V V W - VW, V = V r 



(9) 



This equation gives the integrability condition of (0 and ([8]) 
that must be satisfied by the mass function M(r, v) in order 
for it to be a solution to the collapsing system. Once a function 
M(r,v) solving equation (0 is found the system of Einstein 
equations is fully solved. It is clear that the above equation 
has a non-empty class of solutions to which the self-similar 
collapse models as well as the FLRW class belong. By in- 
tegrating the equation (O we obtain the function t(r,v) that 
describes the time at which the shell labeled by r reaches the 
'epoch' v. Then we can retrieve the singularity curve t s (r), 
namely the curve describing the time at which each shell be- 
comes singular, by setting t s (r) = t(r, 0). 

For all sufficiently regular solutions M(r, v) the singularity 
curve is also regular and can be written near the center as, 



t s {r) = t + Xir + X2r 2 + o(r 3 ) , 



(10) 



where to = t(0, 0) is the time at which the central shell be- 
comes singular and \i = -^rr |t=o- Here xi vanishes due to 
regularity requirements. The tangent to the singularity curve 
close to the center is then determined by the coefficient X2 
and it is then possible to show that it is this coefficient that 
determines the local visibility of the central singularity or oth- 
erwise ( fiUl . HI]). Specifically, if the singularity curve is an 
increasing function of the coordinate r at the center, that is 
X2 > 0, we then have an outgoing null geodesic family com- 
ing out from the central singularity. In that case the singularity 
is locally naked. On the other hand, the singularity is covered 
if X2 <0. 

The FLRW collapse is obtained when we take M = 
mo/v 3 , for which v = v(t) and G becomes ^r- (from which 
we can see that imposing the regularity condition / 2 (0) = 
6mo ensures G = 1 at the center). The singularity curve 
is constant in this case, all shells fall into the singularity at 
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the same time and the event horizon forms before the time of 
formation of the singularity, thus giving rise to a black hole 
final state for collapse. Going to more general case, solving 
the complete second order PDE (O is in general unattainable. 
Nevertheless our purpose here is to extract the crucial infor- 
mation regarding the local visibility or otherwise of the central 
singularity, and this can be achieved by considering a close 
neighborhood of the center. In order to do so we expand Ein- 
stein equations and all relevant functions in powers of r close 
to the origin. 

We therefore studied here a class of spacetimes obtained by 
this procedure and the two-dimensional metric can be given 

as, 

= - f 1 + ) di'W [ 1 + m ' vvi ~ 5my3 r 1 | dr 1 



3m,Q 

(11) 

This line element is a solution of the Einstein equations to 
second order in r which is valid in a small neighborhood of 
r = 0. The metric corresponds to a mass function M(r, v) — 
^§-+m(v)r 2 , a velocity profile f 2 = 6too, and cj)(t(v,r)) 2 = 
Here m(v) is a solution to the corresponding ordinary 
differential equation coming from the expansion of equation 
(|9]i to second order in r, which now becomes 
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(12) 



note that the zeroth order gives the FLRW differential equa- 
tion. 

Once (fT2l is solved we obtain m(v) from which we can 
evaluate the singularity curve as in equation (TTOb to determine 
the coefficient \2, that eventually decides the local visibility 
or otherwise of the singularity. For the above class of solutions 
X2 is given by 



X2 



3 m 



(5r, 



:5m- 
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[ mo 



-dv 
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We note that the 'reality condition' in this case takes the 



form 



+ 



5m t- ' 
3mo 



> and appears at the denominator in \i- 
Therefore those mass profiles m for which it is violated will 
have complex values for \2- Hence we can easily see that 
all possible local behaviours of the central singularity are en- 
tirely determined by the value of \2- The central singularity 
is locally naked if \2 is positive, is covered if it is negative, 
and there is a dispersal of the collapsing central shell if \2 is 
complex. 

Given a certain initial data in the form {to(1), m, v (l)}, a 
solution of equation (TTZb exists and it is possible to solve it 
numerically and obtain the value for \2- If we choose initial 
data sufficiently close to FLRW, meaning with |to(1)| < mo, 
and require m tV (l) to be small to satisfy the 'reality condi- 
tion', we can then evaluate \2 as above and its sign will then 
decide the visibility of the central singularity for the specific 
massless scalar field collapse given by that m(v). 
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FIG. 1. The final outcome of collapse of a massless scalar field de- 
pending on the initial data {m = m(l),m tV — m tV (l)} in a neigh- 
borhood of the initial data for FLRW collapse. For every point in 
the phase space of initial data, the sign of \2 was evaluated numer- 
ically. The thick curve represents the boundary at which \2 — 
that separates the initial data leading to a black hole outcome (on the 
left) from the initial data leading to a naked singularity (on the right). 
Note that m(l) = m, v (l) — corresponds to initial data for the 
FLRW solution. 



What we find then is, for initial data close to FLRW, both 
the collapse outcomes namely the black hole and naked singu- 
larity are equally possible and the space of initial data leading 
to either of them within this specific class is evenly divided 
(see Fig. [T]). The general behaviour found here for the initial 
data sets close to FLRW resembles the results that were found 
in similar scenarios for collapse of perfect fluids and fluids 
sustained only by tangential pressure Q. The main virtue 
of the study of scalar field collapse comes from the fact that 
these are fundamental matter fields obeying the Klein-Gordon 
equation and therefore are possibly valid at all late stages of 
collapse, if quantum gravity effects are ignored. Also, such a 
'stiff fluid' is a fully closed system where there are no global 
free functions, as in the case of perfect fluids or tangential 
pressures when no equations of state are specified. So the 
physical meaning of the matter model comes out in a clear 
and straightforward manner. 

The above result provides intriguing perspective on the 
genericity aspect of occurrence of black holes and naked sin- 
gularity final states in gravitational collapse. The definition 
of genericity in general relativity is a delicate matter that re- 
quires a proper and deeper understanding of both the measure 
and topology of the parent space than is available presently. 
No precisely well-defined criteria are available for the same 
in the gravitation theory today. What we have showed here is 
that there exist classes of solutions for which initial data arbi- 
trarily close to that of the FLRW models lead the collapse to a 
naked singularity. Furthermore, the set of initial data leading 
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to a naked singularity within this particular space of solutions 
does not reduce to a set of zero measure, and the same is true 
for the black hole outcomes. In this sense, within this class, 
the naked singularity formation can be considered somehow 
generic. 

Finally we would like to emphasize that what we deduced 
here is the occurrence of a locally naked singularity that dis- 
proves the stronger version of the censorship conjecture. It is 
possible in principle that the singularity is only locally naked 



and trajectories do come out, but then they eventually all fall 
back into the horizon, without being globally visible, and 
the weaker version of the censorship conjecture may still be 
obeyed. This needs to be examined separately. 
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